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WMAP observations have accurately determined the position of the first two peaks and dips in 
the CMB temperature power spectrum. These encode information on the ratio of the distance to 
the last scattering surface to the sound horizon at decoupling. However pre-recombination processes 
can contaminate this distance information. In order to assess the amplitude of these effects we use 
the WMAP data and evaluate the relative differences of the CMB peaks and dips multipoles. We 
find that the position of the first peak is largely displaced with the respect to the expected position 
of the sound horizon scale at decoupling. In contrast the relative spacings of the higher extrema 
are statistically consistent with those expected from perfect harmonic oscillations. This provides 
evidence for a scale dependent phase shift of the CMB oscillations which is caused by gravitational 
driving forces affecting the propagation of sound waves before recombination. By accounting for 
these effects we have performed a MCMC likelihood analysis of the location of WMAP extrema to 
constrain in combination with recent BAG data a constant dark energy equation of state parameter 
w. For a flat universe we flnd a strong 2a upper limit w < —1.10, and including the HST prior 
we obtain w < —1.14, which are only marginally consistent with limits derived from the supernova 
SNLS sample. On the other hand we infer larger limits for non-fiat cosmologies. From the full 
CMB likelihood analysis we also estimate the values of the shift parameter R and the multipole 
la of the acoustic horizon at decoupling for several cosmologies to test their dependence on model 
assumptions. Although the analysis of the full CMB spectra should be always preferred, using the 
position of the CMB peaks and dips provide a simple and consistent method for combining CMB 
constraints with other datasets. 
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I. INTRODUCTION 

Cosmic Microwave Background (CMB) observations 
have provided crucial insights into the origin and evolu- 
tion of present structures in the universe [l!, [2, |3| . Phys- 
ical processes occurred before, during and after recombi- 
nation have left distinctive signatures on the CMB. The 
most prominent feature is a sequence of peaks and dips 
in the anisotropy power spectrum, the remnant imprints 
of acoustic waves propagating in the primordialphoton- 
baryon plasma at the time of decoupling [J, |a, Q . This 
oscillatory pattern carries specific information on several 
cosmological parameters [Tj . As an example the angular 
scale at which these oscillations are observed provides a 
distance measurement of the last scattering surface to the 
sound horizon at decoupling, hence a clean test of cosmic 
curvature [S]. 

WMAP observations have accurately detected the 
peak structure of the CMB power spectrum. These 
data have constrained the geometry of the universe to 
be nearly flat and have precisely determined other cos- 
mological parameters Q . On the other hand constraints 
on dark energy are less stringent, this is because its late 
time effects leave a weaker imprint of the CMB which is 
diluted by degeneracies with other parameters. Indeed 
other cosmological tests can be more sensitive to the sig- 
nature of dark energy, nonetheless they still require ad- 
ditional information from CMB to break the parameter 
degeneracies. As an example CMB constraints are usu- 



ally combined with those from SN la luminosity distance 
data. Alternatively the CMB can be used in combina- 
tion with measurements of the baryon acoustic oscilla- 
tions (BAO) in the galaxy power spectrum J^. In fact 
the same acoustic signature present in the CMB is also 
imprinted in the large scale distribution of galaxy, thus 
providing a complementary probe of cosmic distances at 
lower redshifts. 

A likelihood analysis of the CMB spectra is certainly 
the more robust approach to implement CMB constraints 
with those from other datasets. This can be very time 
consuming, henceforth one can try to compress the CMB 
information in few measurable and easily computable 
quantities. Recent literature has focused on the use of 
the shift parameter R, and the multipole of the acous- 
tic scale at decoupling la [ll|, ]^ ■ However these quan- 
tities are not directly measured by CMB observations, 
they are inferred as secondary parameters from the cos- 
mological constraints obtained from the full CMB like- 
lihood analysis. Consequently their use as data can po- 
tentially lead to results which suffer of model dependen- 
cies as well as prior parameter assumptions made in the 
analysis from which the values of R {la) have been in- 
ferred in the first place. In contrast the multipole lo- 
cation of the CMB extrema can be directly determined 
from the observed temperature power spectrum through 
model-independent curve fitting. These measurements 
can then be used to constrain cosmological parameters 
provide that pre-recombination corrections are properly 



taken into account. 

In this paper we analyse in detail the cosmological in- 
formation encoded in the position of the CMB extrema 
as measured by WMAP. Our aim is to provide a simple 
and unbiased method for incorporating CMB constraints 
into other datasets which is alternative to that of using 
R and/or la [III, U^ ■ Firstly we estimate the amplitude 
of pre-recombination mechanisms that can displace the 
location of the CMB extrema with the respect to the an- 
gular scale of the sound horizon at decoupling. In partic- 
ular we show that the WMAP location of the first peak is 
strongly affected by such mechanisms, while the displace- 
ments induced on the higher peaks and dips are smaller. 
By accounting for these effects we perform a cosmological 
parameter analysis and infer constraints on dark energy 
under different prior assumptions, including the cosmic 
curvature. We then combine these results with measure- 
ments of BAG from SDSS and 2dF data ^, and con- 
front the inferred constraints with those obtained using 
SN la data from the Supernova Legagy Survey [ij|. Fi- 
nally we test for potential model dependencies of R (and 
la) by performing a full likelihood analysis of the WMAP 
spectra for different sets of cosmological parameters. 

The paper is organized as follows: in Section |TT] we 
review the physics of the CMB acoustic oscillations. In 
Section Hill we discuss the relative shifts of the multipoles 
of the WMAP peaks and dips. In Section [TVl we present 
the results of the cosmological parameter inference using 
the location of the CMB extrema in combination with 
BAO. In Section|V]we confront the results with the SN la 
likelihood analysis from the SNLS sample. We discuss the 
results on the shift parameter in Section IVII and present 
our conclusions in Section IVlIl 



II. CMB ACOUSTIC OSCILLATIONS 

The onset of acoustic waves on the sub-horizon scales 
of the tightly coupled photon-baryon plasma before re- 
combination is natural consequence of photon pressure 
resisting gravitational collapse. The properties of these 
oscillations depends both on the background expansion 
and the evolution of the gravitational potentials associ- 
ated with the perturbations present in the system. In the 
following we will briefly review the basic processes which 
affect the propagation of these waves before decoupling. 
Interested readers will find more detailed discussions in 
[g, |7| . Let consider the photon temperature fluctuation 
00 = AT (monopole), following Hu and Sugiyama 6] its 
evolution is described by 
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where the dot is the derivative with respect to conformal 
time, R = ?>phl^Pf is the baryon-to-photon ratio, k is 
the wavenumber, Cs = c/-y/3(l -I- R) is the sound speed 
of the system with c the speed of light. The source term 
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represents a driving force, where <& and ^ are the gauge- 
invariant metric perturbations respectively. 

It is easy to see from Eq. ^ that the homogeneous 
equation (F = 0) admits oscillating solutions of the form, 



eS°"(7y) = Ai cos krs{r]) + -^ smkr,{T]) 
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where Ai and A2 are set by the initial conditions and 
rs(?7) = Jq Cs{r]')dri' is the sound horizon at time rj. At 
time of decoupling 77* , the positive and negative extrema 
of these oscillations appear as a series of peaks in the 
anisotropy power spectrum. Their location in the multi- 
pole space is a multiple integer of the inverse of the an- 
gle subtended by the sound horizon scale at decoupling, 
namely l^"''' = mla with m = 1,2, ... and 
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where 2;* is the recombination redshift and r{z) the co- 
moving distance to z. 
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with iJg the Hubble constant, |f7/f| = —K/Hq with K 
the constant curvature, f{x) = sin(x),sinh(a;),a; for K > 
0, < and = respectively, and I{z) = J^ dz'Ho/H{z'). 

Scales for which the monopole vanishes also contribute 
to anisotropy power spectrum. In such a case the sig- 
nal comes from the non- vanishing photon velocity Qi 
(dipole) which oscillates with a phase shifted by 7r/2 
with the respect to the monopole Q . Therefore photons 
coming from these regions are responsible for a series of 
troughs in the anisotropy power spectrum at multipoles 
ifP ^nla with n = m + 1/2. 

The full solution to Eq. ([1]) at decoupling reads as [la] : 
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dij'll + R{Tj')f^sin[krs{v*) ~ kr,{ij')]F{ij'), 
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where A3 is set by the initial conditions. As we can see 
from Eq. ([6]) including the driving force F induces a scale 
dependent phase shift of the acoustic oscillations, which 
is primarily caused by the time variation of the gravita- 
tional potential <&. In fact perturbations on scales which 
enter the horizon at the matter-radiation equality expe- 
rience a variation of the expansion rate which causes a 
time evolution of the associated gravitational potentials. 
This mechanism is dominant on the large scales and is re- 
sponsible for the so called early Integrated Sachs- Wolfe 
(ISW) effect [l3|. The overall effect is to displace the 
acoustic oscillations with the respect to the pure har- 
monic series. For a spectrum of adiabatic perturbations 
we may expect this displacement to become negligible on 
higher harmonics since the gravitational potentials decay 
as $ ex (fcry)"^ on scales well inside the horizon. This is 
not the case if active perturbations were present on such 
scales before the epoch of decoupling. 

In order to account for these pre-recombination effects 
a realistic modeling of the multipole position of the CMB 
maxima and minima is given by [l7[ 



Irn = la{m - Ifim), 



where m — 1^2, .. for peaks, and m = 3/2, 5/2, .. for dips; 
(Pm parametrizes the displacement caused by the driving 
force. Because of the scale dependent nature of the driv- 
ing effect discussed above, it is convenient to decompose 
the correction term as ipm = ^ + Sipm, where (p = cpi 
is the overall shift of the first peak with respect to the 
sound horizon, and 6ipm is the shift of the m-th extrema 
relative to the first peak [1§ |. 

It is worth noticing that while the position of the CMB 
extrema depends through la on the geometry and late 
time expansion of the universe, their relative spacing de- 
pends through ip„i only on pre-recombination physics. 



III. PHASE SHIFT OF WMAP PEAKS AND 
DIPS 

WMAP observations have provided an accurate deter- 
mination of the CMB power spectrum. The multipoles 
of the CMB extrema have been inferred using a func- 
tional fit to the uncorrelated band powers as described 
in [1^ . Hinshaw et al. [3| have applied this method to the 
WMAP-3yr data and found the position of the first two 



peaks and dips to be at li = 220.8±0.7, 1 
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412.4±1.9, 



I2 = 530.9 ± 3.8 and /g/s = 675.2 ± 11.1 respectively. 

We want to determine whether these measurements 
provide any evidence for driving effects affecting the 
acoustic oscillations. In order to do so we evaluate the 
relative spacings between the WMAP measured m-th and 
77i'-th extrema. 
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(7) FIG. 1: WMAP spacings of Z3/2, h and I5/2 relative to h 



and the propagated errors cta 



(black solid circles) and propagated errors. The values ex- 
pected from the harmonic series are Ai,3/2 = 1/2, Ai,2 ~ 1 
and Al 5/2 = 3/2 (open circles). Vertical dashed lines delimit 
the expected interval of variation of the relative spacings ob- 
tained by including the shift corrections as parametrized in 
[la | and evaluated over a conservative range of cosmological 
parameter values (see text). The dotted vertical lines include 
the effect of three massless neutrinos. 



Let first consider the spacings relative to the loca- 
tion of the first peak. We find Ai,3/2 = 0.87 ± 0.01, 
Ai,2 = 1.40 ±0.02 and Ai,5/2 = 2.06'±0.05 respectively. 
These estimates are shown in Figure [1] (black solid cir- 
cles), where we also plot the relative spacings as expected 
from a sequence of perfect acoustic oscillations (open cir- 
cles). It is evident that the WMAP inferred values of 
Ai^m lie many sigmas away from those expected from 
the harmonic series. This provides clear evidence that 
the position of the first peak is largely affected by the 
driving force at decoupling. Such a large displacement is 
most likely caused by the early ISW, although an addi- 
tional contribution from isocurvature fluctuations [20| or 
active gravitational potentials [2l| cannot be excluded. 

Let focus now on the displacement of the second peak 
relative to the first one, since Ai_2 > 1 it follows that 
If > S(p2- This implies that the overall shift of h with 
the respect to la is larger than the shift of I2 relative to 
h. As discussed in the previous section this is consistent 
with having the gravitational potentials inside the sound 
horizon scaling as <& oc (fcr;)^^, thus inducing a weaker 
driving force. This can be seen more clearly in Figure [21 
where we plot A3/2,2, A2,5/2 and l^zi^^hji- 

Apart ^2,-ii2 = 0.29 ± 0.01, whose value suggests the 
presence of a non-negligible driving effect still on the scale 
of the first dip, we may notice that all other spacings 
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are statistically consistent with the prediction of the har- 
monic series. 

Therefore these results suggest the existence of a scale 
dependent phase shift of the CMB acoustic oscillations. 
The effect is larger on the scale of the first acoustic peak, 
while it is weaker for the higher harmonics. The upcom- 
ing Planck mission will map more accurately the location 
of the higher peaks and dips and provide a cleaner detec- 
tion of this shift. 

Indeed driving effects are well accounted for by the 
CMB theory as incorporated in standard Boltzmann 
codes 22] . For instance a standard adiabatic spectrum of 
initial density perturbations leads to phase shifts which 
are consistent with those we have inferred here. To 
show this we have used the fitting formulas provided in 
[l_q| for adiabatic models which parametrize ipm in terms 
of the total matter density flmh^, the baryon density 
Q.bh? , the dark energy density at decoupling ^j^^ and the 
scalar spectral index Us- Assuming ^j^% = we evalu- 
ate these formulas over the following range of parameter 
values, 0.08 < n^^h'^ < 0.11, 0.020 < flbh'^ < 0.024, 
0.92 < Us < 1.1 and infer the corresponding intervals for 
the relative spacings A^.m'- These are drawn in Figure[T] 
and [2] as vertical dashed lines. It can be seen that these 
intervals are statistically consistent with the measured 
spacings. Including the contribution of three massless 
neutrinos (dotted vertical hues) shghtly shifts the Ai_m 
intervals further from the expected values of the perfect 
harmonic oscillator. This is because the presence of rela- 
tivistic neutrinos extend the radiation era and therefore 
leads to a more effective early ISW effect on the large 



scales. In contrast we find no differences for the intervals 
of the other peaks and dips spacings. 



IV. PARAMETER INFERENCE 

We perform a Markov Chain Monte Carlo (MCMC) 
likelihood analysis to derive cosmological parameter con- 
straints using the measurements of the WMAP extrema 
discussed in the previous section. Again we account 
for the shift corrections by evaluating the model predic- 
tion for Im using Eg. ([7]), with the displacements ipm 
parametrized as in [l8|. We compute the recombination 
redshift z* using the fitting formulae provided in [23j . 
Cosmological constraints derived from the location of the 
CMB peaks have been presented in previous works (e.g. 
[2J, [23, l2y|). Here our aim is to derive bounds on dark 
energy which are independent of Supernova la data and 
rely only on the cosmic distance information encoded in 
the angular scale of the sound horizon as inferred from 
the multipole position of the WMAP peaks and dips, and 
BAO measurements. 

First we consider flat models with dark energy 
parametrized by a constant equation of state w. We 
then test the stability of the inferred constraints by ex- 
tending the analysis to models with non- vanishing curva- 
ture, flk / 0. We also consider flat dark energy models 
with a time varying equation of state parametrized as 
w = wo + wi{l-a) (CPL) [13, nil. We want to remark 
that for models with wi ^ 1, the dark energy density 
can be non-negligible at early times. Therefore in order 
to consistently account for the shifts induced on the lo- 
cation of the CMB peaks and dips, we compute for each 
model in the chain the corresponding value of ^j^^ so as 
to include its value in the shifts fitting formulae. 

The credible intervals on the parameters of interest 
are inferred after marginalizing over h, flbh'^ and Us 
respectively. We let them vary in the following inter- 
vals: 0.40 < h < 1.00, 0.020 < Qbh^ < 0.024 and 
0.94 < Us < 1.10. Marginalizing over these parameters 
is necessary due to the parameter degeneracies in rx , f's 
and to properly account for the shift corrections fm- 

As complementary dataset we use the cosmic distance 
as inferred from the BAO in the SDSS and 2dF sur- 
veys [13]. These measurements consists of the ratio 
rs(z*)/Z?y(z), where Dv{z) is a distance measure given 

by 



Dv{z)= [{l + zfDA{z)cz/H{z)] 



1/3 



(9) 



with Da{z) — r^f (z)/(l -|- z) the angular diameter dis- 
tance at z. In particular Percival et al. [13[ have found, 
£iy(0.35)/Dy(0.2) = 1.812 ± 0.060. 

In order to reduce the degeneracy with the Hubble pa- 
rameter we also infer constraints assuming a Gaussian 
HST prior h = 0.72 ± 0.08 0. In Figure [3] we plot the 
marginalized 1 and 2a contours in the ^m-w, w-Qk and 
wq-wi respectively. The upper panels correspond to con- 
straints inferred from WMAP extrema alone, while the 




FIG. 3: Marginalized 1 and 2(t likelihood contours from WMAP extrema (upper panels) and in combination with BAO (lower 
panels). Dashed lines correspond to contours inferred under HST prior. The dotted lines in the upper left panel correspond to 
limits inferred assuming Qt = 0.023 and Us — 0.96. 



lower panels include the BAO data. Dashed contours 
are inferred under the HST prior. To be conservative we 
only quote marginalized 2ct limits. We now discuss these 
results in more detail. 



A. Limits from CMB peaks and dips 

As it can be seen in Figure [3] (upper left panel) the 
CMB extrema alone poorly constrain the Qm — w plane. 
In particular the 1 and 2a regions are larger than those 
obtained from the WMAP analysis [^. This is because 
due to the late ISW effect more information about dark 
energy is contained in the full CMB spectrum than just 
in the distance to the last scattering surface as encoded in 
the position of the CMB peaks and dips. Besides several 
degeneracies with other parameters are strongly reduced. 
A direct consequence of this is that our limits on w are 
unbounded from below. After marginalizing over all pa- 
rameters we find r^m = 0.29±Q-23 and w < —0.18 at 2cr. 
A model with J7,„ = 1 is consistent at 95% confidence 
level with the location of the WMAP extrema provided 



that h w 0.42. This is in agreement with the results 
presented in [SO]- On the other hand imposing an HST 
prior (dash contours) reduce the degeneracy in the flm-w 
plane, and the marginalized 2cr limits are flm = 0.16±Qjf 
and w < —0.25 respectively. The upper limit on w im- 
proves if a strong prior on flbh^ and n^. is assumed (dotted 
contours in the upper left panel). As an example impos- 
ing nbh'^ = 0.0223 and n^ = 0.96, we find w < -0.65 at 
2cr. Indeed using the analysis of the full CMB power spec- 
trum provides better constraints. For instance in Fig. [4] 
we plot the 1 and 2cr contours inferred from a MCMC 
likelihood analysis of the WMAP-3yrs spectra in combi- 
nation with the HST prior. The limits are more stringent 
than in the previous case. This is because the amplitude 
of the first peak as well as the relative amplitude of the 
other peaks are particularly sensitive to flm, ^b and h. 
Hence degeneracies contributing to the uncertainties in 
the rim — w plane are further reduced. As mentioned 
before, a robust dark energy parameter inference needs 
the analysis of the full CMB spectrum. However in the 
case one aims to infer constraints from other datasets 
such as SN la or BAO and include CMB information in 




FIG. 4: Marginalized 1 and 2a likelihood contours inferred 
from the full WMAP-3yrs spectra. 



a rapid and simple manner, the position of the CMB ex- 
trema provides a very efficient tool. In fact while the 
CMB power spectrum analysis requires the solution of 
the Boltzmann equation for a given cosmological model, 
the evaluation of the position of the CMB peaks and 
dips only is a semi-analytical computation. As an exam- 
ple running publicly available Boltzamann codes [22| on 
a CPU at 2.3GIIz requires about one minute to compute 
the spectra of a single model, and even using an MCMC 
sampling the overall likelihood analysis still require about 
one hour to reach full convergence of the MCMC chains, 
while using the CMB extrema only takes few minutes. 

In Fig. [3] (central upper panel) we extend our analysis 
of the CMB peaks and dips to non-flat models. Allowing 
for a non- vanishing curvature increases the geometric de- 
generacy and consequently leads to larger uncertainties 
in w. For instance the 2a marginalized constraints are 
w < —0.34 and fix = —0.01 ± 0.05 respectively, and do 
not improve significantly under the HST prior. 

The position of the CMB peaks and dips alone does 
not provide any insight on the time variation of dark en- 
ergy. As it can be seen in Fig. [3] (right upper panel) the 
contours in the wq-wi plane are spread over a large range 
of values. After marginalizing we find wq < —0.55 and 
Wi < 1.68 at 2(7. It is worth mentioning that for increas- 
ing values of wi , dark energy becomes dominant at earlier 
times. In such a case the presence of a non- negligible dark 
energy density at recombination modifies the position of 
the CMB peaks and dips primarily through its effect on 
the size of the sound horizon at decoupling. Therefore 
the location of the CMB extrema (after having accouted 
for the relative shifts) can put an upper bound on the 
time evolution of the equation of state at high rcdshifts 
(i.e. wi). Our analysis shows that in order to be con- 
sistent with the observed peak structure, large positive 
values oiwi ^ 1 are excluded (see also Section IV]). This 
is consistent with the fact that the analysis of the full 
CMB spectrum limits the amount of dark energy density 



at recombination to be less than 10% (otherwise it would 
strongly affect the amplitude and location of the CMB 
Doppler oscillations), hence providing a stringent upper 
bounds on the value of the dark energy equation of state 
at early time (see [3l|, |32[). In contrast models with large 
negative values oi wi < leave no imprint at high red- 
shifts, since in this case the dark energy density rapidly 
decreases for z > 0. Consequently the likelihood remains 
unbounded in this region of the parameter space. 



B. Combined constraints from CMB extrema and 
BAO 

The baryon acoustic oscillations in the galaxy power 
spectrum provide a cosmic distance test at low redshifts. 
Therefore in combination with CMB measurements they 
can significantly reduce the cosmological parameter de- 
generacies. In Fig. [3] (lower left panel) we plot the com- 
bined 1 and 2(7 contours in the flm-w plane. At 95% 
confidence level we find il„i = 0.12 ±0.12 and w < —1.10 
respectively. Imposing the HST prior further constraints 
the dark energy equation of state, w < —1.14. These 
results are compatible with those found in [l3| . A model 
with flrn = 1 is now excluded with high confidence level 
since the combination of CMB extrema and BAO con- 
strain the Hubble parameter in the range h — 0.71 ± 0.20 
at 2(7 (see also [20] )• Interestingly the ACDM case 
(w — —1) appears to be on the edge of the 2a limit, 
hence favoring non-standard dark energy models. Indeed 
unaccounted systematics effects in the BAO data can be 
responsible for such super-negative values of w. On the 
other hand if confirmed this would provide evidence for 
an exotic phantom dark energy component [33| | or in- 
terpreted as the cosmological signature of a dark sector 
interactions (e.g. [SJ])- 

The credible regions for non-flat models are shown 
in Fig. [3] (central lower panel). In this case we find 
f^A' = -0.011 ± 0.064 and w < -0.46 at 2a. These 
bounds do not change significantly under the HST prior. 
In Fig. [3] (lower right panel) we plot the 1 and 2a con- 
tours in the wq-wi plane. Also in this case the bounds 
on a time varying dark energy equation of state remain 
large. For instance we find the marginalized 2a limits 
to be 1^0 < —0.74 and wi < 1.6. Necessarily inferring 
tighter bounds on wi will requires the combination of 
several other datasets such as SN la luminosity distance 
measurements [35|, which is the topic of next Section. 



V. CONSTRAINTS FROM SN lA 

Here we want to compare the results derived in the 
previous Section with limits inferred from luminosity dis- 
tance measurements to SN la. We use the SN dataset 
from the Supernova Legacy Survey (SNLS) [lj|, and for 
simplicity we limit our analysis to flat models. The re- 
sults are summarized in Fig. [5] and Fig. \E\ were we plot 



the 1 and 2a contours in the fim — w and wq — wi planes 
respectively. The shades regions correspond to limits in- 
ferred by combining the SN data with the location of the 
CMB extrema and assuming a hard prior on the baryon 
density and the scalar spectral index, Q.^ — 0.023 and 
Us ~ 0.96 respectively. We have verified that the con- 
straints do not change significantly assuming different 
prior parameter values. 

Let first focus on Fig. [51 We can see that the degen- 
eracy line in the f2m — w plane is almost orthogonal to 
that probed by CMB and BAO, and indeed using the 
SN data requires external information to extract tighter 
constraints on dark energy. A common procedure is to 
assume a Gaussian prior on ^Im consistently with the pa- 
rameter inference from CMB and large scalar structure 
measurements, or alternatively to combine the SN anal- 
ysis with BAO or the CMB shift parameter. Here we 
derive limits by combining the SN data with the position 
of the CMB peaks and dips. This breaks the parameter 
degeneracy, thus providing smaller "credible" contours 
(shaded contours). In particular after marginalizing, we 
find ^m = 0.24 ± 0.11 and w ^ -1.01 ± 0.29 at 2ct 
respectively. We can notice that these limits are only 
marginally consistent with those inferred using BAO in 
the previous Section, thus indicating a potential discrep- 
ancy between the BAO measurements obtained in [1J| 
and the SNLS data [l|. 

Let now consider the case of a time varying equation 
of state. It is obvious that the parameter degeneracy be- 
tween the matter density and the dark energy equation 
of state is increased when additional equation of state 
parameters which accounts for a possible redshift depen- 
dence are included in the data analysis. This can be 
clearly seen in Fig. [S] were we plot the 1 and 2a contours 
in the wq — wi plane. Nevertheless the SN data, differ- 
ently from the case of BAO data in combination with 
CMB extrema (see lower left panel in Fig. [3]), constrain 
Wq in a finite interval. This is because SN la observa- 
tions by testing the luminosity distance over a range of 
redshift where the universe evolves from a matter domi- 
nated expansion to one driven by dark energy, are sensi- 
tive to at least one dark energy parameter (i.e. w or 
Wq) [36|. In such a case adding external information 
breaks the internal degeneracy and leads to finite bounds 
on both dark energy parameters. For instance including 
the position of the CMB peaks and dips, the root-mean- 
square value and standard deviation for wq and wi de- 
rived from the MCMC chains are wq = —1.04 ±0.33 and 
wi = —0.27 ± 2.27 respectively; the best fit model being 
Wq — —1.02 and wi = 0.04. These results are consistent 
with those from other analysis in the literature (see e.g. 
[12]). 



VI. SHIFT PARAMETER 

The geometric degeneracy of the CMB power spectrum 
implies that different cosmological models will have sim- 




FIG. 5: Marginalized 1 and 2a contours in the Qm — w plane 
from SNLS data (solid lines) and in combination with the 
location of the CMB extrema (red and yellow shaded regions). 




FIG. 6; As in Figure [5] in the wq — w\ plane. 



ilar spectra if they have nearly identical matter densities 
flnih^ and flbh^, primordial spectrum of fluctuations and 
shift parameter R = y/flrnHgrKiz^c) [33]- The authors 
of 12] have suggested that since la is nearly uncorrelated 
with R, then both parameters can be used to further com- 
press CMB information and combined with other mea- 
surements in a friendly user manner. For minimal exten- 
sion of the dark energy parameters the inferred values of 
R and la do not significantly differ from those inferred 
assuming the vanilla ACDM model [111, II4I . Indeed dif- 
ferences may arise if additional parameters, such as the 
neutrino mass, the running of the scalar spectral index or 
tensor modes are considered [11]. We extend this analy- 
sis to other models. In particular by running a MCMC 
likelihood analysis of the full WMAP-3yrs spectra we in- 
fer constraints on R and la for models with an extra- 
background of relativistic particles (characterized by the 
number of relativistic species A"e// 7^ 3) [3^, neutrino 
mass |39j , a time varying equation of state parametrized 
in the form of CPL, and a dark energy component with 



Model 


R 


la 


ACDM 


1.707 ±0.025 


302.3 ±1.1 


«;CDM (c^np = 1) 


1.710 ±0.029 


302.3 ±1.1 


wCDM {c^K=0) 


1.711 ±0.025 


302.4 ±1.1 


ACDM nil, >0 


1.769 ± 0.040 


306.7 ±2.1 


ACDM N^ff + 3 


1.714 ±0.025 


304.4 ±2.5 


ACDM Hfc 7^ 


1.714 ±0.024 


302.5 ±1.1 


w{z)QV>M. CPL (c^„^ = 1) 


1.710 ±0.026 


302.5 ±1.1 


ACDM + tensor 


1.670 ±0.036 


302.0 ±1.2 


ACDM + running 


1.742 ±0.032 


302.8 ±1.1 


ACDM + running + tensor 


1.708 ± 0.039 


302.8 ±1.2 


ACDM + features 


1.708 ± 0.028 


302.2 ±1.1 



TABLE I: The 68% C.L. limits on the shift parameter R and 
the acoustic scale derived from the WMAP data. A top-hat 
age prior 10 Gyrs < to < 20 Gyrs is assumed. 



as the DGP scenario [43| deserves a separate comment. 
In these models not only the Hubble law differs from the 
standard ACDM, but also the evolution of the density 
perturbations can be significantly different. Therefore 
unless the evolution of the linear perturbations is under- 
stood well enough as to allow for a precise calculation of 
the CMB and matter power spectra, the use R and la, or 
alternatively of the position of the CMB extrema or the 
distance measurements from BAO might expose to the 
risk of completely wrong results. 



VII. CONCLUSIONS 



perturbations characterized by the sound speed c|,^. We 
also consider models with a running of the scalar spectral 
index, with a non-vanishing tensor contribution (see e.g. 
[40|) and, finally, with extra- features in the primordial 
spectrum due to a sharp step in the inflaton potential as 
in [U. 

As we can see from Table |T] the constraints on R and 
la are stable under minimal modifications of the dark en- 
ergy model parameters, differences are smaller than few 
per cent including the case of a clustered dark energy 
component (c|)^ = 0)- In contrast the confidence inter- 
val of la is shifted by few per cent in the ACDM model 
with the neutrino mass or an extra background of rela- 
tivistic particles, while the values of R are slightly mod- 
ified for a running of the primordial power spectrum or 
the contribution of tensor modes. These results confirm 
previous analysis [ll|, [I4I . 

Although the values of R and la are nearly the same for 
the dark energy models we have considered, this should 
not be considered as an incentive to use these parame- 
ters without caution. For instance there is no specific 
reason as to why one should use the values of R and la 
inferred from the vanilla ACDM, rather than those ob- 
tained accounting for the neutrino mass. Consequently 
one may infer slightly different bounds on the dark energy 
parameters depending whether neutrinos are assumed to 
be massless or not. Moreover the fact that WMAP data 
constrain R and la to be nearly the same for simple dark 
energy models is because the effect of dark energy on the 
epoch of matter-radiation equality and the evolution of 
the density perturbations remains marginal. This might 
not be the case for other models, such as those for which 
the dark energy density is a non-negligible at early times. 
Since this effect is not accounted for in the values of R and 
la inferred from the vanilla ACDM, their use may lead 
to strongly biased results. In contrast the location of the 
CMB extrema is applicable to this class of models as well 
[l8|. A similar consideration applies to inhomogeneous 
models in which the late times dynamics and geometry 
departs from that of the standard FRW universe [42 1. 

The applicability to models of modified gravity, such 



The multipoles of the CMB extrema can be directly 
measured from the WMAP spectra and used to com- 
bine CMB information with other cosmological datasets. 
Corrections to the location of the CMB peaks and dips 
from pre-recombination effects need to be taken into ac- 
count for an unbiased parameter inference. Here we have 
shown that the position of the first peak as measured by 
WMAP-3yrs data is strongly displaced with the respect 
to the actual location of the acoustic horizon at recom- 
bination. This displacement is caused by gravitational 
driving forces affecting the propagation of sound waves 
before recombination. These effects are smaller on higher 
harmonics, indicating the presence of a scale dependent 
phase shift which becomes negligible on scales well inside 
the horizon. 

We have performed a cosmological parameter inference 
using the position of the WMAP peaks and dips in combi- 
nation with recent BAO measurements and derived con- 
straints on a constant dark energy equation of state under 
different model parameter assumptions. 

The method we have presented here is alternative to 
using the shift parameter R and/or the multipole of the 
acoustic horizon at decoupling la- We have tested for 
potential model dependencies of R and la by running a 
full CMB spectra likelihood analysis for different class 
of models. Indeed for simple dark energy models the in- 
ferred constraints on R and la do not differ from those in- 
ferred assuming the vanilla ACDM. Nevertheless we have 
suggested caution in using these secondary parameters as 
data, since hidden assumptions may lead to biased results 
particularly when testing models which greatly depart 
from the ACDM cosmology. 

Indeed we do advocate the use of the full CMB spec- 
tra, particularly for constraining the properties of dark 
energy. In fact more information on dark energy is en- 
coded in the full CMB spectrum than just in the distance 
to the last scattering surface. Nevertheless we think that 
using the location of the CMB extrema provide a fast 
and self-consistent approach for combining in a friendly 
user way the CMB information with complementary cos- 
mological data. 
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